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We develop a phenomenological hydrodynamic theory of coherent magnetic precession coupled 
to electric currents. Exchange interaction between electron spin and collective magnetic texture 
produces two reciprocal effects: spin-transfer torque on the magnetic order parameter and the Berry- 
phase gauge field experienced by the itinerant electrons. The dissipative processes are governed 
by three coefficients: the ohmic resistance, Gilbert damping of the magnetization, and the "/3 
coefficient" describing viscous coupling between magnetic dynamics and electric current, which stems 
from spin mistracking of the magnetic order. We develop general magnetohydrodynamic equations 
and discuss the net dissipation produced by the coupled dynamics. The latter in particular allows 
us to determine a lower bound on the magnetic-texture resistivity. 

PACS numbers: 72.15.Gd,72.25.-b,75.75.+a 



I. INTRODUCTION 

Conduction electrons moving in a ferromagnet interact 
with the magnetization through the exchange interaction. 
If the exchange field is strong and slowly varying in space 
and time, the electron spin will adiabatically follow the 
direction of the magnetization. We may then consider 
electrons with spins up and down along the magnetiza- 
tion direction as two distinct species of particles, and for 
convenience call them spin up/down electrons. As is well 
known, a spin up/down electron wave packet acquires a 
Berry phase^ that influences their orbital motion. In ef- 
fect, the electrons experience a Lorentz force due to "fic- 
titious" electromagnetic fields which are local functions 
of the magnetizationP 

In this fictitious electrodynamics, spin up/down elec- 
trons have opposite charges and different conductivities. 
Their motion and associated currents interact with the 
magnetization through what is commonly called current- 
driven spin-transfer torques. We call this interplay be- 
tween spin currents and magnetization spin magneto- 
hydrodynamics, in analogy to the classical theory of 
magnetohydrodynamics ™ where the magnetic fields cou- 
ple to electric currents in conducting fluids, and the cur- 
rents in turn generate magnetic fields. In our spin magne- 
tohydrodynamics, the Maxwell's equations for the mag- 
netic field are replaced by the Landau-Lifshitz-Gilbcrt 
(LLG) equation for the magnetization. In this paper, we 
neglect full dynamics of the real electromagnetic fields, 
focusing on the spin-related phenomena. 

The electron spin follows the magnetization direction 
perfectly only in the limit of an infinitely large exchange 
field. In reality, there will be some misalignment and 
associated spin relaxation. This is usually described 
phenomenologically as a dissipative spin torq ue w ith a 
coefficient f3 in the Landau-Lifshitz equation! ' 5 ' ' In a 
one-dimensional ring geometry, we will derive the com- 
plete set of coupled spin-magnetohydrodynamical equa- 
tions, starting from the semi-phenomenological dynami- 
cal equations for noncquilibrium currents and magneti- 
zation. We recast the reactive spin torque mediated by 



the Berry phase in this thermodynamic context. In our 
theory, we take an alternative view that the (3 term arises 
from a correction to the Berry-phase electromotive force 
(EMF) in the equation of motion for the charge current, 
with the appropriate dissipative spin torque established 
by the Onsager reciprocity. 

This physics is presently vigorously studied (exper- 
imentally as well as theoretically) in the contexts of 
curren t-driven m agnetic excitations and domain-wall 
motiorP^ 7 ' 8 ' 9 ' 10 ! and the reciprocal spin accumula- 
tions and voltages generated by the fictitious gauge 
fields IIlI 12 l 1 3 | i4 | i5 lT6] g mce t ne mesoscopic regime (mainly 
dealing with variants of magnetic spin valves, tunnel 
junctions, and magnetic multilayers) is at present well 
explored,^ we will limit our attention here to the case of 
continuous magnetic systems. 



II. NONDISSIPATIVE SPIN TORQUE 

Since the underlying physics is rich and complex in 
the most general setting, we will limit our discussion to 
a simple setting, which we believe captures all the es- 
sential ingredients of the spin magnetohydrodynamics. 
Consider a uniform current in a ferromagnetic ring, as- 
suming for simplicity incompressible electric flows (the 
continuity equation prohibits current inhomogeneities for 
an incompressible electron fluid). The electric current is 
then the only dynamical variable describing the electron 
fluid. The magnetic texture here could be a domain wall 
or magnetic spiral, for example (in higher dimensions we 
could have topological twists and kinks such as vortices, 
hedgehogs, or skyrmions). See Fig. [T for a schematic 
of the setup. In the Landau-Lifshitz phenomenology of 
ferromagnetic dynamics well below the Curie tempera- 
ture, only the instantaneous direction of the magnetiza- 
tion m(x,t) (or, equivalcntly, spin density) is assumed 
to be a dynamic variable. The magnitude of the spin 
density S along m is assumed to be uniform and con- 
stant in time. We will separately drive the current with 
a time-dependent external magnetic flux <f>(t) inside the 
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FIG. 1: (color online). Schematics of our principal "study 
case:" Uniform electric current J(t) carried by itinerant elec- 
trons can be driven by the external magnetic flux <&(t) gen- 
erating the EMF £ = — 9 t $/c. The magnetic texture m(x,t) 
responds to the effective field H(a;, t), which may have an ex- 
ternal contribution applied to the wire independently of 
The reactive magnetohydrodynamic coupling stems from the 
Berry phase which is acquired by the electron spin (shown 
in blue) following the instantaneous magnetic profile (shown 
in red) around the loop. $' corresponds geometrically to the 
solid angle enclosed by the electron spin. Coupled dissipative 
processes arise once we relax the projection approximation, 
allowing for some orientational spin mistracking and dephas- 
ing as electrons propagate through the magnetic texture. 



ring, and the magnetic dynamics with a magnetic field 
h(:r, t) applied directly to the wire. 

The first step in our phenomenology is to identify 
the free energy T as a function of the thermodynamic 
variables J and m(x, t) (or their thermodynamic con- 
jugates), which completely determine the macroscopic 
state of our system, assuming local thermal equilibrium. 
Neglecting spin, the gauge-invariant free energy associ- 
ated with an electric current in the ring is given by 
F(J,$) = (J -$/c) 2 /2L, where we define LJ to be 
the current corresponding to the canonical momentum 
of the electrons. L is the self-inductance of the ring and 
c is the speed of light. However, spin up/down electrons 
propagating through a quasistatic magnetic texture^ ac- 
cumulate also a Berry phased which gives a fictitious 
contribution to the vector potential associated with a 
fictitious EMF P] This vector potential is given (in some 
convenient gauge) A' x = (hc/e)sin 2 (8/2)d x (j), pro- 

ducing gauge-invariant fictitious flux, 



= I dxA' x 



he 



<p dx(l — cos 9)d x 4> . (1) 
2e J 



{9, <fi) are the spherical angles parametrizing m(x). e > 
is minus the electron charge. Eq. ([I I is the flux associated 
with spin-up electrons adiabatically following magnetic 
texture, with the opposite result for spin-down electrons. 

The free energy accounting for the Berry phase be- 
comes 

P(J, &[m{x, i)]) = [J-(® + p&)/cf /2L , (2) 



where p is the polarization of the spin s-dependent con- 
ductivity a s : p = (cr-f — a[)/(a^ + a{) (assuming fast 
spin relaxation or halfmctallic fcrromagncts) . The elec- 
tric current is given by 

J = -ccV = [J- ($ +p&)/c] /L = djT' , (3) 

which is thus the thermodynamic conjugate of J . The 
equation of motion for current in our simple electric cir- 
cuit is given by Ohm's law, 



d t J = Ld t J + d t ($ + p&)/c=-RJ. 



(4) 



where R is the resistance of the wire. Naturally, the 
dynamic Berry phase is seen to give a contribution to 
theEMFP 



£' = -pdt&'/c ■ 



V <j> dxm ■ (d x m x dtm) . 



(5) 



which is a well-known resultP (We defined V = ph/2e.) 

Now that the free energy of the current is coupled to 
the magnetization of the ring through the Berry-phase 
flux, there will be a corresponding reactive coupling of 
the magnetization to the current. We describe magnetic 
dynamics by the Landau-Lifshitz-Gilbcrt equation 19 



dtm = H x m/S — am x <9 t m , 



(6) 



where the effective field H is defined by the functional 
derivative, H = d m !F (so that locally H _L m), and 



a is the dimensionless Gilbert damping^ parameter. 
The total free energy of our magnetoelectric system is 
J"(m, J, $) = _F(m)+JF'( l 7, $'[m(jc, £)]), where T(m) 
is a standard free energy of the ferromagnet. Variation of 
the T' with respect to m gives current-driven spin torque 
applied to the magnetic dynamics)^ r' = d m T' x m, 
where d m T' = d m ^'d^T' = -pJd m $'/c. Differentiat- 
ing Berry phase ([!]) with respect to m, we find 

t' = VJd x m . (7) 

Since H/2e is the electron spin-charge conversion factor, 
we can give another interpretation of this term, ft is sim- 
ply the rate of change of the angular momentum of the 
conducting electrons with spins locked to the magnetic 
profile. The spins of the up/down electrons rotate in the 
opposite directions so that, if the spin up/down conduc- 
tivities are the same (and thus V = 0), the net change in 
their angular momentum vanishes. Putting this term on 
the left-hand side, we get 

d t m — VJd x m/S = ^^(m) x m/S — am x d t m . (8) 

The left-hand side of this equation is the rate of change 
of the total angular-momentum density of the magneto- 
electric systempl while the right-hand side gives the usual 
LLG torque on the system. 
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III. DISSIPATIVE SPIN TORQUE 

LLG equation ^ with torque ^ and Ohm's law Q 
with the fictitious EMF §5§ now constitute coupled equa- 
tions of our spin magnetohydrodynamic theory with the 
reactive coupling mediated by Berry phase ([TJ. We re- 
produce them here for clarity (after putting the magne- 
tization equation in the Landau-Lifshitz form): 



d t J 



-RJ , d t m 



H x m aH 

(l + a 2 )S 



(9) 



These are the equations of motion for a quasistationary, 
thermodynamic system near equilibrium. 22 In equilib- 
rium, the current J is zero and magnetization is static. 
Out of equilibrium, the first-order time derivatives of 
(J , m) are completely specified by the instantaneous val- 
ues of their thermodynamic conjugates (J, H). The right- 
hand side is a linear expansion in these conjugates with 
dissipative coefficients R and a that cause the system to 
relax back to equilibrium. So far, the dissipation in the 
current and magnetization is separate and physically un- 
related. We now add the dissipative couplings which will 
be key results of this paper. 

We proceed phenomenologically by adding to the cur- 
rent equation Q correction AS' to the Berry-phase EMF 
and correction R' to resistance, due to coupling with the 
magnetic texture m(x,t). The modified Ohm's law then 
becomes: 



BtJ = -(R + R')J + AS' 



(10) 



To avoid a slew of uninteresting coefficients and 
anisotropics, we will constrain the phenomenology by as- 
suming spin-rotational symmetry of the magnetic texture 
and the inversion symmetry of the wire. Under the lat- 
ter, m — > m, J — ► —J, d x — ► —d x , and £' — » —£'. In the 
spirit of the standard quasistationary description J22l we 
expand only up to the linear order in the nonequilibrium 
quantities J and dtm, so that terms of the form, e.g., 
J 2 dtm - d x m are excluded. To the second order in d x m, 
the only possible terms satisfying these requirements are: 

AS' - R'J = 

/(3 2 V 2 f 
dx d x m ■ d t m — n — — J i> dx(d x m) . (11) 
aS J 

The first term stems physically from a spin mistracking 
of electrons propagating through the magnetic texture™ 
Since the mistracking should scale as 1 / A xc (vanishing in 
the limit of infinite exchange A xc ), we may anticipate the 
dissipative coupling to be governed by a small parame- 
ter (3 ~ H/t s A xc , where t s is a characteristic (transverse) 
spin-dephasing time. The r\ term in Eq. ( 10 I describes the 



resistance associated with magnetic texture, which is of- 
ten discussed in the context of magnetic domain walls.^ 
Both terms in Eq. (Ill are odd under time reversal, like 



that including in Eq. (Ill a reactive term of the form 
([5]) would not add anything new to the following consid- 
erations, as long as we treat V as a phenomenological 
coefficient. 

Our modification of Ohm's law must respect the 
Onsager reciprocity principle. 22 Substituting dtm from 
Eqs. (|9| into Eq. (JllJ, we see how the effective field H 
(which is conjugate to m) affects the dynamics of J . The 
Onsager theorem is now readily applied to determine how 
the electric current J (which is conjugate to J) should 
modify the dynamics of m. We write the final result as 
a correction to the spin torque ([7]): 



At' = 13V Jm x d x m . 



(12) 



The complete equation of motion of the magnetic texture 
in the LLG form thus becomes 

d t m= H x m/S-am x d t m + At'/S, (13) 

with t' implicitly included in H. 



Eqs. (10) and (13 1 are our final coupled deterministic 



equations. We can rewrite them in a more explicit form 
as 

Ld t J + (R + R')J + d t <S>/c = 

V j> dx d x m ■ {(3 — mx)<9 f m , 
5(1 + amx)<9 t m + m x H = VJ(l + flmx)d x m . (14) 

Here, the deterministic spin-torque contribution ([7]) is for 
clarity separated out of the effective field H, which here 
consists of the usual purely magnetic contributions. The 
left-hand sides in these equations contain the ordinary 
Ohm's law (corrected for the magnetic-texture resistance 
R') and the LLG terms, respectively, while the right-hand 
sides describe the reactive Berry-phase coupling and its 
dissipative f3 correction. 



ohmic resistance and Gilbert damping. Finally, we note 



Eq. ( 12 1 was derived microscopically in Rcfs. 4 6 24 25 , 
relating [3 to electron spin dephasing: (3 ~ h/T s A xc (con- 
sistent with our anticipation above). Its Onsager coun- 
terpart in Eq. (Ill was first obtained phenomenologi- 
cally in Ref. [TJ] and microscopically in Ref. [T3] These 
11 (3 terms" are now accepted to be crucial in understand- 
ing current-driven magnetic dynamics and the reciprocal 
gauge fields. 



IV. DISSIPATION POWER 

Suppose we perturb our system with some nonequi- 
librium current and magnetic texture, after which the 
system evolves back toward equilibrium according to the 
equations of motion, producing entropy. If the system is 
steadily driven, the heat will be dissipated to the envi- 
ronment at some finite rate. From standard thermody- 
namics, the dissipation power is 
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P[m(x,t),J(t)] = -Jd t J ~ i dxll ■ d t m = RJ 



dx 



B 2 T> 2 1 
aS(d t mf -2(3VJd x m-d t m + r]^-—- J 2 (d x m) 2 . (15) 

ab 



r 



According to the second law of thermodynamics, the dis- 
sipation (151 must always be positive, which means that 
?y > 1. This gives us the lower bound on the resistivity 
of the magnetic texture: 

P = V — 5-(3s m , 
ab 



R 2 V 2 
2 > ^V(5,m) 2 
ab 



(16) 



In models where a comes solely from the coupling of 
the magnetization to the conducting electrons (which is 
in fact believed to be the dominant cause for Gilbert 
damping in metallic ferromagnets) , we may expect the 
lower bound (16 1 to give an estimate for the texture re- 
sistivity. For a mean-field Stoner-model treatment of 
Gilbert damping, we found a = f3, while for an s — d 
model we had a — (s/S)f3, where s is the portion of 
spin density carried by the s electrons, S is the total spin 
density, and (3 = H/t s A xc in both cases (with the spin- 
dephasing time r s governed by the magnetic and spin- 
orbit impurities) P In both models, therefore, ab = s(3, 
giving for the resistivity estimate (up to the second order 
in spatial derivative) 



P>{f3V 2 /s){d x mf 



(17) 



which involves only quantities related to conducting elec- 
trons. Taking parameters relevant to Permalloy wiresf^ 
p ~ 1, (3 ~ 10~ 2 , domain-wall width of 20 nm, and 
the mag netization of 10 3 emu/cm 3 , we find the resistiv- 
ity (17) to be p ~ 1CP 4 /ifi-cm. This is smaller than 
the domain-wall resistivity calculated to the (1/A XC ) 2 
order in spin mistracking of the magnetic profile (but 
still quadratic order in texture), in the absence of spin 



relaxation p2l whose overall prefactor appears to be larger 
than in our Eq. (17 1 for transition metals. We thus con- 
clude that our r\ may in practice be much larger than 
unity (which is the lower bound necessary for the consis- 
tency of our phenomenology) . 

Let us also note in the passing that in the special case 
of a = (3 and 77 = 1, the magnetic dissipation (151 ac- 
quires a very simple form: 



P[m(x,t)] 



ab <j> dx 



VJ_ 
b 



ctm 



(18) 



which is nothing but the Gilbert dissipation with the ad- 
vective time derivative D t = dt + vd x (v = —VJ/S). 
It is clear that this limit describes dissipative magnetic 
dynamics that are simply carried by the electric flow at 
speed v. In this case, the spin torques disappear if we 
write the LLG equation ^ with D t in the place of dt . 8 
V. THERMAL NOISE 

At finite temperatures, thermal agitation causes fluctu- 
ations of the current and magnetization, which are cor- 
related due to their coupling. A complete description 
requires that we supplement the stochastic equations of 
motion with the correlators of these fluctuations. It is 
convenient to regard these fluctuations as being due to 
a stochastic external magnetic field Sh and a stochastic 
current source 8 J: their noise correlators are then related 
to the dissipative coefficients of the theory according to 
the fluctuation-dissipation theorem (FDT). Constructing 
the noise sources by following the standard procedure 
our final coupled stochastic equations become: 



Ld t J + R(J+5J) + d t <S>/c = V j dxd x m ■ (/3 - mx)9 ( m , 
5(l + amx)a t m + mx (H + Sh) = V Jd x m + V(J + 6 J) (3m X d x m, 



(19) 
(20) 



where we have explicitly separated the deterministic spin-torque contribution VJd x m out of the effective field H, 
which here consists of the usual purely magnetic contributions. The left-hand sides in these equations contain the 
ordinary Ohm's law (corrected for the magnetic-texture resistance: R = R + R') and the LLG terms, respectively, 
while the right-hand sides describe the reactive Berry-phase coupling and its dissipative (3 correction. 
Writing {J, H} = —7 ® {d t J, d t m}, we read out the "matrix" 7 from Eqs. ( 19 1 and (20 1: 



1 

i? 7 



—d x m, 7h(x),j= -^4m, 



(3 2 V 2 

lhi(x) : h z ,(x<) = be tl 3 mj(x)6(x - x) + aS5 U '5(x- x) —d x mi(x)d x mi, (x) 



(21) 



where t % ^ k is the antisymmetric Levi-Civita tensor. Symmetrizing matrix 7 immediately produces Langevin sources 
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satisfying the FDT, ^ in the limit that fiw < k B T: 



{6J{t)6J{t')) = 2k B TS(t - t') /R , (6J{t)Sh(t')) = , 
{8h i {x)8h v {x')) = 2k B T \aS5 u >5(x - x') - ((3 2 V 2 / 1 Rjd^d^m,, S(t - t') 



(22) 



Apart from the obvious contributions, we have a magnetic field noise proportional to (3 2 , in the form of a nonlocal 
tensor Gilbert damping. The nonlocal Gilbert damping is apparent, if the electrons are not externally driven, dt<& = 0, 
in the limit L — > of a large ring, in which case the magnetic equation decouples to give 

V 2 / 

5(1 + amx)9 t m + m x (H + Sh + Sh') = + /3mx)<3 x m <J> dx' d x >m ■ (/? — mx)9 f m . (23) 

Here, we moved the spin torque driven by the Nyquist noise to the left as 

Sh! = -VSJm x d x m . (24) 

Sh' thus enters the equation as a statistically independent current-driven noise source. Writing the right-hand side of 
Eq. d23l as 



m x <p dx'~K(x, x')d t m(x') . 



where 



V 



K iV {x,x') = — (m x d x m - (3d x m) i (m x 9 x /m + pd^m)^ , 



and extracting the symmetric part of the tensor Kai(x,x'), we arrive at the total Gilbert damping tensor 
G u >(x,x') = aSa'Six - x') + — ^ [(m x d x m)i(m x d^m)^ - (3 d x mid x imi-] . 



(25) 



(26) 



(27) 



r 



This is exactly the form required by the FDT, consistent 
with the correlator for Sh + Sh.' . The effective Gilbert 
damping can thus appear both negative and positive in 
different regions. The minimal texture resistivity (16 1, 
however, insures that we have a nonnegative damping 
globally. This Gilbert damping originates physically in 
the spin torques that are generated by the magnetically- 
driven fictitious EMF. Nonlocal d x d x i magnetic noise was 
recently constructed in Ref. [55] (neglecting spin relax- 
ation and P) by heuristically converting Nyquist current 
noise into magnetic fluctuations via adiabatic spin trans- 
fer. Although the DFT-required nonlocal d x d x i Gilbert 



tensor (27) was established in that paper (apart from the 
1 piece), only here we are able to derive it directly from 
the fundamental Langevin sources of the coupled mag- 
netohydrodynamic theory, dictated by the FDT. As es- 
timated in Ref. 26, this nonlocal contribution to Gilbert 
damping is in practice important (in comparison to a) in 
nanoscale magnetic structures. 



VI. SUMMARY 

We developed a general phenomenological theory of 
magnetohydrodynamic coupling in isotropic metallic fer- 



romagnets. The reactive coupling between magnetic tex- 
ture dynamics on the one hand and electric flows on the 
other stems from the Berry phase accumulated by elec- 
tron spin following the quasistationary magnetic texture. 
Dissipative terms of the coupled dynamic equations orig- 
inate in the electron spin mistracking of the magnetic or- 
der parameter and the associated spin dephasing. Apart 
from the usual Gilbert damping, the latter leads to a 
viscous coupling between electric currents and magnetic 
texture dynamics, parametrized by a single parameter (3. 
We also obtain a small correction to the texture resistiv- 
ity at order 1 . Finally, our thermodynamic description 
of the magnetohydrodynamic coupling allows us to de- 
rive the stochastic Langevin contributions to the effective 
field and electric current, according to the fluctuation- 
dissipation theorem. 
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